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Causality creates an asymmetry between space and time, even though the wave equation treats 
them on equal footing. In this work, we leverage this asymmetry to construct a cross-mapping 
between space and time. This cross-mapping is applied to simplify scattering in space-varying 
media, by eliminating the infinite reflections between interfaces. The method is shown to transform 
the implicit transfer matrix method into an explicit method for the analysis of electromagnetic field 
scattering by a stratified medium. 
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I. INTRODUCTION 

Time-varying media alter the angular frequency spec¬ 
trum of waves propagating across them. Such trans¬ 
formations are rarely discussed in electrodynamics text¬ 
books, except for the Doppler effect, which is essentially 
a change in frequency due a moving source or scattering 
medium as seen by a stationary observer. However, time- 
varying media systems not involving any motion, such as 
magnetoplasma and magneto-elastic media, may be more 
practical in electromagnetic applications [l|-y. Related 
effects include frequency conversion, time reversal, and 
angular spectrum engineering possibilities, with a wide 
range of potential applications in electromagnetics and 
optics. 

In this work, we present a space-time cross-mapping 
approach where scattering in space-varying media is 
treated as scattering in corresponding time-varying me¬ 
dia. This approach may simplify the analysis of space- 
varying scattering problems by leveraging the breaking, 
due to causality, of the inherent space-time symmetry 
in the wave equation. Indeed, while the wave equa¬ 
tion admits solutions that correspond to waves propa¬ 
gating forward and backward in space and in time, how¬ 
ever, causality prevents temporal backward-propagating 
waves. In the case of discontinuous time-varying media, 
such abruptly switched media, this breaking of space- 
time symmetry prohibits multiple reflection between the 
temporal interfaces. Treatment of wave scattering in such 
media is inherently simpler than in corresponding space- 
varying media with multiple reflections between the inter¬ 
faces. However, certain challenges also arise in the cross- 
mapped configuration, such as the presence of boundaries 
outside the line-of-sight (LoS) of the source and observer. 
This challenge arises precisely due to space-time asymme¬ 
try, where boundaries outside the LoS in a time-varying 
configuration are not permitted to reflect waves back¬ 
ward in time. A second closely related challenge is the 
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proper cross-mapping of sources from spatial into tem¬ 
poral configurations. 

This work is organized as follows: In Section im the 
principles of wave scattering at temporal discontinu¬ 
ities, which are the simplest form of time-varying me¬ 
dia, are detailed. The relevant continuity conditions, 
conservation laws and an ongoing controversy pertain¬ 
ing to the electromagnetic wave momentum are discussed 
and canonical scattering at simple temporal discontinu¬ 
ities are analyzed. In Section imi the space-time cross¬ 
mapping method is formulated. First, the appropriate 
mapping between the incident and scattered waves in 
temporal and spatial configurations, and the relevant 
mapping of media properties are established. Next, the 
treatment of corresponding boundary and continuity con¬ 
ditions is presented, followed by the mapping of sources. 
In Section lEl applying the cross-mapping is shown to 
simplify the transfer matrix method (TMM) by trans¬ 
forming it from an implicit into an explicit method. Con¬ 
clusions are given in Section O 


II. WAVE SCATTERING AT TEMPORAL 
DISCONTINUITIES 

Wave scattering at spatial discontinuities has been 
thoroughly covered in the literature. Examples of canon¬ 
ical scattering at simple spatial discontinuities include 
the problems of two half-spaces, slabs, and stratified me¬ 
dia 0- Understanding canonical scattering at simplified 
discontinuities provides a basis for understanding more 
involved scattering phenomena. By analogy, understand¬ 
ing wave scattering at temporal discontinuities is the ba¬ 
sis for understanding more involved temporal scattering 
phenomena and may be fundamental for complete under¬ 
standing of electrodynamic phenomena in general. 

In all cases, understanding wave scattering at an inter¬ 
face requires establishing continuity conditions through 
the interface. In what follows, we show the difference be¬ 
tween the continuity conditions at spatial and temporal 
interfaces, and apply these conditions to scattering by 
simple temporal discontinuities. 
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A. Continuity Conditions and Conservation Laws 

In the scattering of electromagnetic fields at a spatial 
interface, the conventional continuity conditions imply 
the continuity of the tangential electric and magnetic 
fields at the interface. Even though the mathematical 
formulation of these conditions is not rigorous in the case 
of discontinuities , it yields satisfactory results away 
from the discontinuity in most practical configurations. 
However, a more rigorous interpretation of the continu¬ 
ity conditions in terms of actual physical principles is still 
preferable, especially when investigating unconventional 
configurations. The physical principles involved are con¬ 
servation laws, such as the conservation of energy and 
momentum across the discontinuity. 

Scattering at passive and lossless spatial discontinu¬ 
ities implies the conservation of energy flow, whence all 
incident power is transformed into scattered power. In 
the case of harmonic waves, the energy conservation con¬ 
dition reads 

(S‘) = (S'-) + (S*), 

where (S) is the time-averaged power density vector, and 
the superscripts i, r and t respectively refer to the inci¬ 
dent, reflected and transmitted waves. In addition to 
the conservation of energy, in linear media the rhythm 
of the wave, i.e. its temporal dependence, is also con¬ 
served to satisfy the phase matching at the discontinu¬ 
ity. This phase matching condition is what dictates the 
change in the wave vector of the scattered wave. Let a 
harmonic wave propagating in the ^-direction be given 
by ex.p{i[±k^z ± cjjt]), where j G {i, r,t}. Assuming an 
incident wave propagating in the positive ^-direction, 
ex.-p{i[k[Z — cjit]), the phase matching condition at a dis¬ 
continuity located at z = 0 reads 

Vt, (1) 

where 7 and r are the coefficients of the scattered waves 
in the incidence medium and scattering medium, respec¬ 
tively. Condition o thus implies the conservation of the 
angular frequency 

LUi=UJr=LUt. (2) 

For simplicity, the signs of ujj are chosen to be the same. 
Consequently, the wavevectors of the reflected and trans¬ 
mitted waves have to change to satisfy the wave disper¬ 
sion relation, = k^c^^ where c is the speed of light in 
free-space. This implies that 

^ = = ( 3 ) 

rii fir rit 

where nj is the refractive index in the corresponding 
medium. 

Alternatively, at a temporal discontinuity, the phase 
matching condition results in conservation of the pattern 
of the wave, i.e. the wavevector in the case of harmonic 


waves. Assuming an incident wave propagating in the 
positive 2)-direction, the phase matching condition at a 
temporal discontinuity occurring at t = 0 reads 

^ ^ (4) 

where ^ and ^ are the coefficients of the scattered waves 
in the medium after the discontinuity. Here, it is essential 
to notice that, due to causality, no wave is permitted to 
reflect in time and propagate into the incidence medium, 
in contrast to the case of Eq. o. The condition (| 1 ]) thus 
implies the conservation of the wavevector 

ki kf k\). ( 5 ) 

Similarly, for simplicity, the signs of k^ are chosen to be 
the same. Consequently, the angular frequencies of the 
scattered waves have to change to satisfy the wave dis¬ 
persion relation. This implies that 

ujiUi = UfUf = —c^b^b- (6) 

It is also worth noting that negative sign in front of cJb 
means that the wave is propagating in the negative 2)- 
direction, since the wavenumber has not changed accord¬ 
ing to Eq. Accordingly, at a temporal discontinuity, 
an incident wave splits into a forward- and a backward- 
propagating waves while incurring shift in the angular 
frequency. 

The conservation of k suggests that corresponding 
physical conservation law invoked at a temporal disconti¬ 
nuity is the conservation of momentum flow of the wave. 
This implies that the momentum flow of the incident 
wave must be equal to the momentum flow of the scat¬ 
tered wave at the discontinuity, which splits the scattered 
wave into a forward- and backward-propagating waves. 
The continuity condition for a space-harmonic wave may 
thus be expressed as 

(&) = (Gb) + (Gf), 

where (G) is the space-averaged momentum density. 

The definition of the momentum density of electromag¬ 
netic waves in a dielectric medium is a controversial is¬ 
sue. Two expressions exist with claimed theoretical and 
experimental support and still no conclusive resolution. 
The two expressions are derived from the photon momen¬ 
tum in a dielectric medium with refractive index n. In 
Ref. 0 , Minkowski derived the (canonical) photon mo¬ 
mentum as Pm = nhwjc^ where h is the reduced Planck 
constant, yielding the electromagnetic momentum den¬ 
sity Gm = D X B, where D and B are the displacement 
field and the magnetic flux, respectively. On the other 
hand, in Ref. [H, Abraham derived the (kinetic) photon 
momentum as pa = hwjinc)^ yielding the electromag¬ 
netic momentum density Ga = x H, where E and 
H are the electric and magnetic fields, respectively. This 
difference in electromagnetic momentum expressions re¬ 
sulted in what is known as the Abraham-Minkowski con¬ 
troversy [i,0- 
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Since the Abraham-Minkowski controversy is still not 
resolved in a definite fashion, the proposed cross-mapping 
method is presented for scalar waves. This scalar rep¬ 
resentation avoids ambiguity in the correct form of the 
continuity condition at temporal discontinuities, since 
the expressions for the energy and momentum densities 
are unambiguous for scalar waves. This representation 
however does not limit the applicability of the proposed 
method, since what is actually modeled is a spatial scat¬ 
tering problem, not a temporal scattering one, and the 
correct continuity conditions are therefore well-defined. 
Hence, the results of the scalar method are directly ap¬ 
plicable to the vector electromagnetic case with proper 
translation as detailed in Section uni 

B. Simple Discontinuities 


CJ 



To understand the space-time cross-mapping method, 
it is necessary to first study scattering at simple tem¬ 
poral discontinuities, such as single interfaces, temporal 
slabs, and stratified temporal media. In what follows, 
we continue with the example of plane-wave incidence on 
a temporal discontinuity in the previous section. The 
phase matching condition implied that k is conserved 
and uj is changed in the scattering medium. To deter¬ 
mine the scattered wave coefficients, the continuity of 
the time derivative of the wave across the discontinuity 
is invoked, which translates to the continuity of the mo¬ 
mentum density in the vectorial case. The continuity of 
time derivative of the wave reads 


- V^. (7) 

Conditions and o with the relations and m 
determine the scattered wave coefficients as 


^ _ ni + rif 
^ “ 2ni 

rij - nf 

^ “ 2ni 


(8) 

(9) 


Figure [T] depicts the frequency shift at simple spatial 
and temporal discontinuities. The slopes of the world 
lines (the straight lines passing through the origin) in the 
uj-k diagram designate the phase velocities in the two 
media. Conservation of uj at the spatial discontinuity, 
as indicated in Eq. ([2j), is designated by the horizontal 
dashed arrows. The k values at the intersection points of 
the arrows with the world lines give the wavenumbers of 
the scattered waves, as indicated in Eq. ([3]). Conservation 
of k at the temporal discontinuity, as indicated in Eq. ®, 
is designated by the vertical arrows. The uo values at 
the intersection points of the arrows with the world lines 
give the angular frequencies of the scattered waves, as 
indicated in Eq. ©• 

A natural extension to the temporal discontinuity is 
the temporal slab. A temporal slab is constructed by 
switching an unbounded medium and then switching this 


FIG. 1. Wave frequency shift at an interface between two 
media. At a spatial interface, the wave number of the inci¬ 
dent wave is shifted (dash-dotted arrow) in the reflected and 
transmitted waves, while at a temporal interface, the angular 
frequency of the incident wave is shifted (dash-double-dotted 
arrow) in the forward- and backward-propagating waves. 


medium back to the initial medium after a time interval 
r in between. Wave scattering through such slab follows 
as a direct extension of the discontinuity case. The waves 
scattered by the second slab interface are causal. There¬ 
fore, they are not allowed to reflect back in time and 
interact with the first interface. This is a manifestation 
of the asymmetry between space and time due to causal¬ 
ity. Assuming a slab of refractive index n 2 and duration 
r, embedded in a background medium with refractive in¬ 
dex ni, the composite (global) coefficients of the forward- 
and backward-propagating waves after passing through 
the temporal slab are found by applying Eqs. (|8]) and (|9]) 
at each one of the two interfaces of the slab, yielding 

5 = [6,16,26-“^" + > (10) 

s = , (11) 

where the subscripts p, q designate the transition across 
a temporal discontinuity from a medium with into a 
medium with nq. Figure [2al shows a pictorial diagram of 
wave scattering through a temporal slab indicating local 
and composite scattered wave coefficients. 

Figure [3] compares wave scattering through spatial and 
temporal slabs. In the case of spatial slabs, the composite 
reflection and transmission coefficients, F and T, respec¬ 
tively, combine an infinite series of reflections and trans¬ 
missions, while in the case of temporal slabs, the com¬ 
posite forward- and backward-propagation coefficients, E 
and E, respectively, are each a sum of only two waves. 
The difference between the two material slabs is due to 
causality, which prevents the waves from reflecting from 
temporal interfaces into the past. It is also worthwhile 
noting that the direction of the trajectories of harmonic 
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waves do not pertain to energy transfer, but rather to 
the direction of phase evolution, as required by the phase 
matching condition (j4]). 

Extending the temporal slab results to stratified (mul¬ 
tilayer) temporal media and temporal photonic crys¬ 
tals is straightforward. The coefficients of forward- and 
backward-propagating waves for each additional tempo¬ 
ral interface are recursively constructed from those at the 
previous one in the same fashion Eqs. m and m are 
constructed from Eqs. m and as shown in Eig. [20 
Thus, the forward- and backward-propagating wave co¬ 
efficients in layer j + 1 are constructed from the those in 
layer j as 

2i+i = , (12) 

5i+i = . (13) 







P - 
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FIG. 2. Wave scattering in temporal layered media. |(a)| In 
a temporal slab, the composite scattering coefficients may be 
directly deduced from the local scattering coefficients at the 
two temporal interfaces of the slab. |(b)| In a general multilayer 
temporal medium, the composite forward- and backward- 
propagating wave coefficients at a temporal interface may be 
directly deduced from the local scattering coefficients and the 
wave coefficients in the prior layer. 


III. CROSS-MAPPING FORMULATION 

As demonstrated in Section [III wave scattering at tem¬ 
poral discontinuities is simpler than that at spatial dis¬ 


continuities due to causality, which prevents reflection 
into the past. The proposed space-time cross-mapping 
method aims at leveraging this asymmetry between spa¬ 
tial and temporal discontinuities to simplify the compu¬ 
tation of scattered waves in a space-varying environment. 
The proposed method thus cross-maps space and time 
coordinates, allowing for the treatment of scattering at 
spatial discontinuities as scattering at temporal discon¬ 
tinuities. That is, the spatial configuration is mapped 
into the corresponding temporal configuration, solved in 
this simpler configuration, and then mapped back to the 
initial configuration. This mapping aims at reducing 
the computational complexity, since no infinite reflec¬ 
tions and transmissions are present in a temporal scatter¬ 
ing configuration with a finite number of discontinuities. 
Moreover, resembling marching-on time schemes 0[i2, 
albeit in space, the mapping may be employed in devel¬ 
oping marching-on in space computational techniques. 

Two particular issues require careful consideration in 
applying the method. The first is accounting for conti¬ 
nuity and boundary conditions at discontinuities outside 
the source-observer LoS, and the second is accounting 
for the presence of sources and their emitted backward- 
propagating waves. These two aspects have no direct 
counterparts in temporal configurations due to causality. 
The boundaries outside the LoS in a temporal configura¬ 
tion have no effect on the scattered wave, since waves can¬ 
not scatter back into the past to reach the observer. Sim¬ 
ilarly, in the case of a source, the backward-propagating 
wave emitted by the source has no direct equivalent in a 
temporal configuration. 

In the following sub-sections, first, space-time cross¬ 
mapping is formulated for source-free configurations 
without discontinuities outside the LoS. Next, the cross¬ 
mapping method is extended to account for the presence 
of discontinuities outside the LoS and the presence of 
sources, which completes the formulation of the method. 


A. Space-Time Cross-Mapping 

Consider the two cases of plane-wave scattering by a 
spatial material slab, shown in Eig. [30 and by a tem¬ 
poral material slab, shown in Eig. [30 In the spatial 
configuration, the incident wave propagates in the pos¬ 
itive 2 )-direction, while, in the temporal configuration, 
the incident wave is considered to propagate in the neg¬ 
ative ^-direction. The geometry of the two scattering 
configurations are seen to be identical if one considers 
the composite wave coefficients rather than the individ¬ 
ual trajectories, and inverts the direction of all the wave 
trajectories. This suggests that a mathematical mapping 
may by established between the two configurations. 

The space-time cross-mapping procedure is thus car¬ 
ried out, with the aid of Eig. [0 as follows: first, the 
spatial configuration coordinate axes are mapped as 2 : ^ 
—cf and ct z'. Second, the wave trajectories are re¬ 
versed, resulting in mapping the transmitted wave in the 
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FIG. 3. Dielectric slab cross-mappin g be tween spatial config¬ 
uration and temporal configuration. |(a)| Spatial configuration 
of the dielectric slab showing the incident, the infinite number 
of reflected and transmitted waves, a nd th eir composite repre¬ 
sentations, /, r, and T, respectively ! (b)| Temporal conhgura- 
tion of the dielectric slab showing the incident, two forward- 
and two- backward-propagating waves, and their comp osite 
representations, S, and S, respectively. Note that in |(a)[ 
the inci dent wave is propagating in the positive ^-direction, 
while in |(b)| the incident wave is propagating in the negative 
; 2 ;^-direction. The arrows in both illustrations pertain to the 
evolution of the phase. 


spatial configuration into the incident wave in the cross- 
mapped temporal configuration. Third, the scattered 
waves in the cross-mapped temporal configuration are 
determined employing temporal treatment. Fourth, the 
substitution n ^ 1/n in applied to the final expressions 
of the composite wave coefficients. This substitution is 
necessary to account for the difference in frequency shift 
between spatial and temporal configurations, as shown in 
Eqs. (|3)) and (|6j). Finally, the obtained composite coeffi¬ 
cients of the scattered waves are normalized with respect 
to composite forward-wave coefficient, S. Table [I] pro¬ 
vides a summary of cross-mapping for source-free config¬ 
urations. 

To illustrate the cross-mapping procedure, consider the 
case of scattering through a dielectric slab of width I with 
refractive index n 2 embedded in a background medium 


with refractive index ni. Following Table [H and using 
Eqs. ([TQ|) and ([TT]) lead to the following relation 


\ _ ^2g2in2/cZ 
^ [1 ~ 

_ ^ _ g2in2/cZj 


1 

T’ 

F 

T’ 


where 7 = (ni — n 2 )/{ni + 77 , 2 ). Now F and T may be 
directly computed without needing to sum up an infi¬ 
nite number of terms. The spatial scattering process has 
thus been simplified, by leveraging the asymmetry due 
to causality. 

Note that space-time cross-mapping is not limited to 
harmonic waves and may be extended to pulsed waves 
in a straightforward fashion. A pulsed wave may be 
expanded in terms of its temporal spectral components 
through Fourier transform, assuming all media are linear, 
then the cross-mapping is carried out for each spectral 
component of the pulse. Finally, spectral composition 
through inverse Fourier transform of the spectral compo¬ 
nents of the scattered waves is carried out to determine 
the composite scattered waves. 

As previously mentioned, cross-mapping is not re¬ 
stricted to scalar waves. It may directly be translated 
to vectorial electromagnetic configurations in the cases 
of longitudinally-invariant propagating modes and trans¬ 
verse electric-magnetic (TEM) fields. The cross-mapping 
in the vectorial case requires the substitution n ^ Z 
instead of n 1/n, where Z is the impedance of the 
medium in the case of TEM fields, or the transverse 
impedance in the case of longitudinally invariant modes. 


B. Discontinuities outside the Line-of-Sight 

Discontinuities existing between the incident wave and 
the observer are naturally accounted for after applying 
the cross-mapping. However, if a discontinuity exists 
outside the LoS connecting the incident wave and the 
observer, it will not be directly accounted for in the cross- 
mapped temporal configuration. The reason is causality, 
where a discontinuity outside the LoS is not permitted 


TABLE 1. Cross¬ 
figurations 

Spatial 

z 

ct 

n 

I 

n 

TI 

F 

T 


mapping between spatial and temporal con- 


Temporal 

1/n 

-I'E 

-f 

T/h 
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to reflect waves back in time towards the observer. Spe¬ 
cial consideration must therefore be given to such cross¬ 
mapping configurations. 

Consider the example of the dielectric slab of the previ¬ 
ous section, where now the observer is located inside the 
slab at distance d away from the exit interface, as shown 
in Fig. 03 Following the cross-mapping procedure out¬ 
lined in the previous section and simply cross-mapping 
the trajectories from the incident wave to the observer 
yields here an incorrect result, since the second inter¬ 
face of the slab, formed in the future of the observer, 
is not accounted for. The proper cross-mapping has to 
be carried out between proper boundaries, namely inter¬ 
faces where the values of the wave are prescribed. In this 
case, the boundaries are radiation boundaries located at 
^ ^ Too. This requirement is consistent with the na¬ 
ture of harmonic waves as they exist, by definition, in all 
space. The composite coefficients of the scattered wave 
at the observer are found by normalizing the scattered 
wave coefficients to E following Table [I] as 


p{d) = ^ + 

^ 2n2 


p{d) =e^ri2kd2^ 


in2kd^‘^d _ in2kd^‘^ 

2n2 


T, 


where the appropriate refractive index substitution has 
been performed, and the phase factors pertain to the lo¬ 
cation of the observer inside the slab. 

Continuing with the observer inside the slab example, 
consider the case shown in Fig. [40 where the slab is 
backed by a hard boundary at its exit interface. The hard 
boundary is equivalent to a perfect electric conductor 
(PEC) boundary in the case of electromagnetic waves. 
Following the requirement that cross-mapping has to be 
carried out between boundaries, proper cross-mapping 
treats the hard boundary as a source that emits forward- 
and backward-propagating waves that exactly annihilate 
each other at the hard boundary. The cross-mapping 
procedure is otherwise carried out as in the previous case. 
The composite coefficients of the forward- and backward- 
propagating waves inside the slab are found as 


p{d) = 


p{d) = 


^—in2kd 


^in2kd 

^^ 2e^ri2kl _^^_^^^^in2kl ' 


C. Sources 

Waves emitted by a source always propagate away from 
the source, which in the one-dimensional configuration 
under consideration are the forward and backward di¬ 
rections. A discontinuity existing outside the LoS con¬ 
necting the source and the observer will scatter the wave 
emitted by the source back towards the observer. In 



FIG. 4. Spa tial slab configuration with cross-mapped wave 
trajectories. |(a)| Proper trajectory mapping starts from the 
radiation point at z ^ oo if an open boundary exists, (b) 


Proper trajectory mapping starts from a hard boundary be¬ 
having as a source. 


a cross-mapped temporal configuration, such scattered 
wave is not permitted to reflect back in time and hence 
will not reach the observer yielding incorrect evaluation 
of the waves at the observer. 

Simple cross-mapping of waves emitted by a source is 
not possible due to causality, since the waves emitted by 
the source and propagate away from the observer, ac¬ 
cording to Sec. IIII Ai should be mapped to waves that 
propagate backward in time. Proper source modeling 
thus implies wave splitting at the source with an apriori 
undetermined source strength a. The cross-mapping is 
carried out as previously described until the final bound¬ 
ary, whether it is a hard boundary or a radiation bound¬ 
ary. Applying the proper boundary condition at the final 
boundary determines the cross-mapped source strength 
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a. The corresponding source condition in the cross- 
mapped configuration is thus given by 


[Giz-,zXz:,-=0, 

(14) 

§-JG{z-,z')]X':=a, 

(15) 


where Eq. 0 establishes the continuity of the Green’s 
function, G, across the source, and Eq. establishes 
the jump of the derivative of G across the source. In a 
spatial configuration, a = 1, i.e. the source strength is 
always unity. 

An example involving a source in the cross-mapped 
configuration is shown in Eig.O where Eig.[5a|is the spa¬ 
tial configuration consisting in a source located at a dis¬ 
tance d away from a dielectric half-space with refractive 
index n, and Eig. [5b] is its corresponding cross-mapped 
configuration. In the cross-mapped configuration, the 
incident wave is split into two waves at the source, an ac¬ 
tual wave and an auxiliary wave, with coefficients a and 
1 — (a, respectively. The auxiliary wave is an anti-causal 
wave that is necessary for proper modeling, but has no 
physical significance otherwise. This auxiliary wave does 
not exist in the spatial configuration, hence its name. To 
determine a, the radiation condition in the dielectric half¬ 
space implies the vanishing of all backward-propagating 
waves in the dielectric as 

+ (1 - a) = 0, 


which yields 


1 

a =-=-. 

2 _ L^‘^ikd 

The wave coefficient observed at z = d + / away from the 
half-space is thus 

-I _ ^ 2ikd 

r J ■ 

Combining the proper source, continuity conditions 
and boundary conditions outside the LoS completes 
the proper formulation of the space-time cross-mapping 
method. 



IV. APPLICATION IN TRANSFER MATRIX 
METHOD 

The transfer matrix method (TMM) is a powerful com¬ 
putational method used to analyze the propagation of 
electromagnetic waves through a stratified medium [a 
However, when there are multiple interfaces, the 
multiple reflections between these interfaces may inter¬ 
fere destructively or constructively. The composite re¬ 
flection or transmission through a multilayer structure is 
the sum of an infinite number of reflections, which scales 
up the computational burden. 



(a) 



(b) 

FIG. 5. A spatial configuration wi th a source located at a 
distance d away from a half-s pace . |(a)| Wave trajectories in 
the spatial configuration, and |(b)| wave trajectories in cross- 
mapped configuration. In the cross-mapped configuration, 
the source splits in the incident wave into two waves, with an 
apriori undefined strength a. The boundary condition at the 
end of the tracing process determines a. The exit boundary 
condition in this configuration is the radiation condition in 
the dielectric at ct' ^ 00 . 


The TMM takes advantage of the longitudinal invari¬ 
ance of each layer of the stratified medium with eigen¬ 
mode expansion in the transverse cross-section, if neces¬ 
sary. Such expansion allows for analyzing inhomogeneous 
structures as long as they can be approximated by a 
stack of layers which may be inhomogeneous in the cross- 
section, while invariant in the longitudinal direction, as 
shown in Eig. ini In such layers, the field at the end of 
the layer may be simply deduced from the field at the 
beginning of the layer through a simple phase accumula¬ 
tion. Additionally, the field continuity condition across 
the boundary between two layers may be expressed as a 
simple matrix operation. Scattering through a stack of 
layers may thus be represented as a matrix systeni, which 
is the product of the individual layer matrices |15l-[l^. 
The computational burden increases with the increase in 
number of layers, hence employing a space-time cross¬ 
mapping may reduce this burden by increasing the com- 
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putational speed or reducing the memory requirement. 



FIG. 6. Schematic of a layered medium or device. The inci¬ 
dent field is designated by /, while the composite reflection 
and transmission coefficients are designated by F and T, re¬ 
spectively. Notice that the cross-section of any given layer 
may be inhomogeneous as long as the layer is longitudinally 
invariant. 


A. Transfer Matrix Method 


Application of the TMM requires each layer of medium 
to be longitudinally invariant. This invariance allows for 
the expansion of electromagnetic field in terms of plane- 
waves, the eigenmodes of the two-dimensional Cartesian 
plane, in the cross-section of each layer. This eigenmode 
expansion may be mathematically represented by dx 
ikx and dy iky^ assuming the ^-direction to be the 
longitudinal direction. Next, the transverse components 
of the wave vector are normalized as = k-Jk^^ with 
j G {x^y}. The longitudinal distance is also normalized 
as z = koz. Substituting these elements into Maxwell’s 
equations and eliminating the longitudinal fields yields 


dz 


^ = 0, 


with the transverse fields vector 


( 16 ) 


and the matrix 


^ = 


Ex 

Ey 

Hx 

^Hy, 


EL = 


0 

0 

kx ky 


0 

0 

Cr -^ 

. Pr 

_ kx ky 

fir 


kx ky 
Cr 

0 

0 


_ kx ky 

Cr 

0 

0 


where and are respectively the effective relative 
permittivity and permeability of the layer considered. 

Equation is a matrix differential equation, whose 
solution takes the form 

(z) = We^^c, 

where c = W“^^(0), ^(0) is the field at the input in- 
ter face of the layer, and W are there eigenvectors and A 
are eigenvalues of EL. 

The continuity conditions at the interface between two 
layers implies the continuity of the tangential fields as 

^i+i (0) = {koLj), 

which may be expressed in terms of the eigenmodes of 
the layers as 

Cj-^i = Tj • cj, (1"^) 

where Tj = exp(Aj/coI/j) is the matrix that 

transfers the field from the layer j into the layer j + 1, 
and Lj is the longitudinal length of the j layer. 

The system matrix is finally constructed from the 
transfer matrices of all layers as 

ct = T • c/, (18) 

where c/ and Ct are the fields at the entry and the exit 
ports of the stratified medium, respectively, and T = 
Tjv • Ttv-i • ... Ti • ^W/ is the composite transfer 

matrix. Equation m describes an implicit scheme, since 
the unknown composite coefficients exist in c/ and ct- 
Solving the matrix equation ([T8]) yields the composite 
reflection and transmission coefficients of the stratified 
medium. 


B. Space-Time Cross-Mapping in TMM 

Space-time cross-mapping may be advantageously in¬ 
tegrated into the TMM with minimal modification and 
with the benefit of reducing the computation burden. 
The resulting method is a marching-on in space method, 
where the scattered field is deduced from the incident 
field without having to account for the infinite reflections 
between the layer interfaces. 

To integrate the cross-mapping method into the TMM, 
the procedure outlined in Section IIII Al is followed and 
c/ and Ct are cross-mapped to Ct and c/, respectively. 
Next, the continuity condition (HZD is replaced by the 
matrix generalization of Eqs. m and m as 

Cj = Tj • Cj+i, ( 19 ) 

where Tj = Wj“^Wj+i exp(Aj+i/coI/j+i) is the cross- 
mapped transfer matrix from layer j+l to layer j, Wj are 
the eigenvectors corresponding to the cross-mapped con¬ 
figuration. The system matrix equation is constructed 
from Eq. m as 


Ct = T • c/. 


( 20 ) 
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where T is the cross-mapped system transfer matrix. The 
marching-on in space nature of the cross-mapped config¬ 
uration yields the implicit scheme in Eq. (|2Q|) . where Ct 
is directly computed from c/, which contains all the com¬ 
posite scattered field coefficients. This is in contrast to 
Eq. m, where the scattered field coefficients exist both 
in c/ and Ct, and the composite scattered field coeffi¬ 
cients are only available after solving the matrix system. 
Since the scattered field is obtained without the need to 
solve a linear matrix system in the cross-mapped TMM, 
the computational burden may be significantly reduced. 

We should also note that the space-time cross-mapping 
method may be also applied to other techniques that 
analyze longitudinally invariant layers, such as rigorous 
coupled-wave analysis (RCWA) and plane-wave ex¬ 
pansion method (PWEM) M- These techniques are in¬ 
herently similar to the TMM, but they do not require 
homogenization of the layers of the medium. 

V. CONCLUSION 

In this work, we have introduced the novel concept of 
space-time cross-mapping with application to wave scat¬ 
tering. The proposed cross-mapping allows for treat¬ 


ing scattering in a space-varying configuration as scat¬ 
tering in a time-varying environment. Cross-mapping 
yields simplified expressions for the composite coefficients 
of the scattered wave by leveraging the asymmetry be¬ 
tween space and time due to causality. Methods to prop¬ 
erly account for boundaries outside the line-of-sight be¬ 
tween the source and the observer, and to properly ac¬ 
count for the presence of sources are presented. Space- 
time cross-mapping is not limited to harmonic waves 
and may be extended to pulsed wave through temporal 
spectrum composition. In addition, the presented cross¬ 
mapping for scalar waves may be directly extended to 
transverse electric-magnetic electromagnetic waves and 
longitudinally-invariant modes. 

The proposed cross-mapping may have potential ap¬ 
plication to existing analytic and computational scatter¬ 
ing schemes, such as the transfer matrix method, rig¬ 
orous coupled wave analysis, and plane-wave expansion 
method. Cross-mapping may be also leveraged to decou¬ 
ple the scattering matrices in coupled filter design. 

Euture work may include dispersion analysis of the 
cross-mapped configuration for computational schemes, 
and extension of the method to three-dimensional con¬ 
figurations. 
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